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Abstract 

We study the effective hydrodynamics of neutral black branes enclosed in a finite 
cylindrical cavity with Dirichlet boundary conditions. We focus on how the Gregory- 
Laflamme instability changes as we vary the cavity radius R. Fixing the metric at 
the cavity wall increases the rigidity of the black brane by hindering gradients of 
the redshift on the wall. In the effective fluid, this is reflected in the growth of the 
squared speed of sound. As a consequence, when the cavity is smaller than a critical 
radius the black brane becomes dynamically stable. The correlation with the change 
in thermodynamic stability is transparent in our approach. We compute the bulk 
and shear viscosities of the black brane and find that they do not run with R. We 
find mean-field theory critical exponents near the critical point. 



1 Introduction and conclusions 



The complex non-linear dynamics of black holes and black branes, governed by Ein- 
stein's equations, can in some regimes be efficiently captured by an effective theory 
for collective degrees of freedom. Some of these degrees of freedom are worldvolume 
embedding coordinates, associated to the elastic dynamics characteristic of brane- 
like objects, but there are also hydrodynamic degrees of freedom — e.g., pressure 
and velocity — associated to horizon dynamics. The hydrodynamic approach dates 
back to [T] and has been greatly extended following the fluid/gravity correspondence 
of [2J. In the context of asymptotically flat black branes, it features within the ef- 
fective blackfold theory of [3j HJ [5] . A framework for relating all these approaches is 
presented in [6]. 

Here we are interested in using the hydrodynamic methods for investigating 
a peculiar effect of black branes: the classical instability of their horizons to the 
formation of ripples along their worldvolume, discovered by Gregory and Laflamme 
in [7]. Refs. (H [8] have shown how this phenomenon is very neatly captured in 
the hydrodynamic theory of blackfolds: it is simply an instability of fluctuations 
of the pressure of the effective black brane fluid, i.e., a sound-mode instability. In 
this paper we investigate it further by introducing an additional parameter to gain 
control over the stability of the system. To do so, we place the black brane inside a 
finite cylindrical cavity of fixed radius R. Then we find the solution for a fluctuating 
black brane, with regularity conditions at the horizon, in the hydrodynamic limit. 
We analyze how the effective black brane theory, in particular its stability under 
hydrodynamic fluctuations, changes as the cavity radius R is varied. This allows us 
to study several issues: 

Correlated instabilities — ghosts vs. tachyons. One reason to expect that 
enclosing the black brane in a cavity should have an effect on its stability comes 
from the Correlated Stability Conjecture (CSC) [9], which links classical dynamical 
stability to local thermodynamical stability. Since it is known that the specific heat 
of the black brane in a cavity changes from negative for cavity radii greater than a 
critical value, R > R c , to positive when R < R c [10], the validity of the conjecture 
requires that the Gregory-Laflamme instability of the black brane disappears at the 
critical radius R = R c . 

To put our study in the right context, it is worth discussing the status and proper 
interpretation of the CSC. This has often been taken as the statement that 
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• (CSC:) translationally invariant horizons have a tachyonic perturbation mode 
if and only if they are locally thermodynamically unstable. 

Indeed, a large part of the studies of the CSC have focused on the presence or 
absence of a static, zero-mode perturbation of finite wavelength (a tachyon) that 
would mark the onset of the instability [HI [12] • In particular, ref. p3] studied 
the static zero-mode for the black brane in a cavity and showed that it disappears 
precisely when the cavity size reaches R = R c . 

This form of the CSC, however, is incorrect: examples of black branes are known 
which are thermodynamically stable but nevertheless have tachyonic instabilities 
[15J. Instead, the local thermodynamical stability of black branes is more appro- 
priately related to the presence of massless ghost excitations, rather than tachyons. 
These two kinds of unstable modes are quite different. For excitations that in some 
range of wavenumbers k have a dispersion relation approximately of the type 

to 2 = c 2 k 2 + m 2 , (1.1) 

we say we have a tachyon when m? < 0. The static zero-mode corresponds to 
k = ko = a/ '—m 2 /c 2 and u — 0. Instead, we say we have a ghost when c 2 < 0, and 
in particular a massless ghost has H = Imu = \/—c 2 k. 

The argument why thermodynamical instabilities of a translation-invariant hori- 
zon are connected to massless ghosts is simple [HI [15], H] . A horizon that is trans- 
lationally invariant can support perturbations of arbitrarily long wavelength. In 
the cases where the frequency of these perturbations vanishes as the wavelength 
diverges, they are hydrodynamic modes, which are either fluctuations of conserved 
quantities or Goldstone modes. Both of them feature in the local thermodynam- 
ics of the fluid. In the case of main interest to us here, the conserved quantity is 
the energy, and it is an old result (which we reproduce in sec. [4] below) that small 
fluctuations in the energy density propagate along the fluid with squared velocity 



where s is the entropy density of the fluid and Cy its specific heat at fixed volume. 
Obviously, a local thermodynamic instability, with Cy < 0, results in unstable 
perturbations of wavenumber k whose amplitude grows exponentially in time like 
exp(^y —v 2 kt). More generally, any local thermodynamic instability of the black 
brane gives rise to a long-wavelength, hydrodynamic instability. In the terms used 
above, we have a massless ghost with c 2 = v 2 < 0. 
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Then, the CSC as stated above must be replaced by a statement of Correlated 
Hydrodynamic Stability: 

• (CHS:) translationally invariant horizons have massless ghost excitations if 
and only if they are locally thermodynamically unstable. The ghost is a long- 
wavelength, low imaginary frequency, hydrodynamic instability of the horizon. 

Since horizons are stable to fluctuations of very short wavelength^ the ghost 
instability at small k must disappear at some larger k = k Q > 0, i.e., lmu(k Q ) = 0. 
If also R,eu(k ) = 0, then this is a zero-mode. In other words, a hydrodynamic 
ghost instability (and hence a local thermodynamic instability) of the horizon will 
typically be accompanied by a tachyonic zero-mode at finite k. This is indeed the 
case for the GL instability of neutral black branes (also when in a cavity). But the 
converse need not be true: a tachyonic instability need not turn into a hydrodynamic 
ghost instability at very long wavelengths^ and hence need not be related to a local 
thermodynamic instability. 

In this paper we construct explicitly the ghost, hydrodynamic unstable pertur- 
bation of the black brane in a cavity, and show that it turns into an oscillatory 
(damped) sound wave when R < R c . 

Increasing rigidity. Our analysis of the effective hydrodynamic theory also gives 
a concrete intuitive picture of why and how the instability disappears as the cavity 
size is reduced. The squared speed of sound of the effective theory at finite R is 
a monotonic increasing function of decreasing R. A larger speed of sound reflects 
a higher rigidity of the system. This comes about because fixing the metric on 
the cavity wall at finite R makes it harder for the geometry to fluctuate, and in 
particular prevents the creation of worldvolume gradients of the redshift on the 
wall. In the hydrodynamic theory, these gradients have the effect of an acceleration 
of the fluid that opposes the creation of inhomogeneities along the worldvolume. 
This works to make the system more rigid and therefore less unstable, until the 
instability disappears. 

Viscosities do not run with R. The solution for the spacetime metric for a 
fluctuating black brane in a cavity of radius R, to first-derivative order in the fluc- 
tuations, allows us to compute the stress-energy tensor of the effective fluid including 

^n very short scales the horizon is indistinguishable from Minkowski space, which in any 
healthy gravitational theory is stable at sufficiently short wavelengths. 

2 For instance, it can become a homogeneous tachyonic mode with Imw(fc = 0) = \J —m? . 
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dissipative effects. The values we obtain for the shear and bulk viscosities (and for 
their ratios to the entropy density) are the same at all values of R. This is very 
likely related to a similar result obtained in the context of black branes in AdSs in 
[16], and may be a feature of a larger class of black branes. 

Spectrum and criticality. The inclusion of dissipative terms in the fluid equa- 
tions gives us an improved approximation for the spectrum of unstable modes. The 
dispersion curves show clearly that the instability weakens as the critical point is 
approached. Although the hydrodynamical theory cannot capture all the physics of 
the critical state, it nevertheless indicates that critical exponents are of mean-field 
theory type, a result which is borne out by the numerical computations of [13J. 

In the remainder of the paper we elaborate on all these points in detail. Sec. [2] 
introduces the black brane in a cavity as a static system. The solution for its 
fluctuations to first-derivative order is discussed in sec. [3j This solution forms the 
basis for the study in sec. [4] of stability from a hydrodynamical perspective, and its 
connection to local thermodynamic stability. Sec. [5] computes the effective viscosities 
of the black brane in the cavity of radius R, and then uses them for obtaining the 
dispersion relation for unstable modes. We conclude with a brief discussion of the 
critical point and the appearance of mean-field critical exponents. 



2 Static black brane in a cylindrical cavity 
2.1 Geometry 

We write the metric of a black p-brane inD = 3 + p + n spacetime dimensions in 
the form 

dv 2 

als 2 = {-f{r)u a u b + P ab ) da a da b + — + r 2 dn [n+1) (2.1) 

f{ r ) 

with 

/(O = 1 - % (2-2) 

and where 

Pab = Vab + u a u b (2.3) 

is the projector onto spatial direcions orthogonal to the timelike vector u a with 
normalization r\ ab u a u h — — 1. 

We put the black brane inside a cylindrical cavity bounded by a 'wall' that 
extends along the brane worldvolume directions a a and which, in the transverse 
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directions, is a sphere S n+1 at finite radius r = R. We denote quantities measured 
on the cavity wall with a caret. The metric induced on the wall is 



h tiu dx fl dx u = h ab da a da b + R l dVt(, n+1 ) 



with 



Kb = -f{R)u a u b + P ab 
= -u a u b + P ab . 



(2.4) 



(2.5) 



Indices of hatted tensors will be raised and lowered with this metric. The velocity 

u a = ^/J{R)u a (2.6) 

is unit-normalized with respect to this metric by absorbing the redshift factor on 
the wall. The orthogonal projector, instead, is not modified, 



ab lab 



(2.7) 



since in the geometry (2.1) the spatial worldvolume directions do not suffer any 



gravitational deformation. 

The geometry of the wall of the cavity is characterized by giving, in addition to 
the induced metric, the extrinsic curvature tensor 



Q,u = -\^KR)d R h 



(2.8) 



Out of this we obtain the Brown- York quasilocal stress-energy tensor on the wall. 
Since we are only interested in the dynamics of the worldvolume, we only consider 
the components of the tensor along the directions a a , and we integrate them over 
the transverse S n+1 of radius R. The result is that 



T, 



n 



ab 



8ttG 
8ttG 



n+1 R n+1 [e ab -h ab e 



-(n + l)R n ^J(R)u a u b + d R (R n+1 ^/J(R~j) P ab 



(2.9) 



2.2 Physical magnitudes 



Eq. (2.9) is a perfect-fluid stress-energy tensor 



Tab = £ U a Ub + PP ( 



ab 



(2.10) 
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with energy density and pressure 

e = --^(n + l)B"V7{R), (2.H) 
P = _ e+ ^ ir » aaym= _ e+ ^^ (, 12) 

In addition, we can assign an entropy density and temperature to the system 

Aixr Q ^f(R) 

The temperature is modified relative to its asymptotic value by the redshift factor at 
the wall, but the entropy density does not depend on R: it is obtained as s = S/V, 
where the total entropy S is computed from the horizon area, and the spatial volume 
V does not undergo any variation as R changes. 

The system satisfies the thermodynamic Euler relation 

e + P = fs, (2.15) 

and the first law 

de = fds, (2.16) 

for variations that keep fixed the cavity radius R. 

For the record, we note that when R is allowed to vary, the first law becomes 

de = Tds — <j w da w (2-17) 

where the wall area-density a w and tension a w are 

„= ^^±m. (2.i8) 

When R can vary the wall is regarded as a dynamical object, and one gets the 
coupled dynamics of the black brane/wall system. However, although this might 
be of interest, for the remainder of the paper we will regard the wall only as a 
non-dynamical boundary condition. 

2.3 No subtraction required 

When R — > oo both e and P diverge, owing to the non-compactness of the space. A 
simple remedy to this is to subtract the stress-energy tensor associated to a surface 
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in Minkowski space with the same induced metric h^. However, we do not need 
this for our purposes. The reason is not merely that we keep R finite and thus 
divergences are absent. More important, the intrinsic worldvolume dynamics that 
we are interested in is not affected by the subtraction. A surface at constant r = R in 
Minkowski spacetime has O a ;, = and the stress-energy tensor comes entirely 
from the curvature of the S n+1 of radius R. Then 

T^ = ^(n + l)R n h ab . (2.19) 

Since we keep R fixed, this stress-energy tensor is of 'vacuum-type', i. e., proportional 
to the worldvolume metric h ab and with constant energy density. This is inert: it 
lacks any hydrodynamic behavior, which is associated with a breakdown of local 
Lorentz invariance and the presence of inhomogeneities on the worldvolume. 

Therefore, the subtraction does not affect the hydrodynamics of the brane, and 
we shall not implement itj^] 



3 Fluctuating black brane 

We promote the parameters u a and tq in the solution to worldvolume collective 
degrees of freedom, i.e., slowly- varying functions of a a . The remaining parameter, 



R, is kept fixed. Following [2], to the now fluctuating metric (2.1) we add correcting 
functions f^ u such that the total metric 



ds 2 = ( r? a6 H —u a (a)u b (a) J da a da b H + r 2 dfl( n+1) 



+f tlu dx^dx u (3.1) 

is a solution to the field equations. We fix the radial coordinate by choosing it to 
be orthogonal to the worldvolume and normalized to measure the S' n+1 -area-radius. 
Then Jq^ = 0. Working to leading order in derivatives, the correcting functions 

f^dx^dx" = f ab da a da b + 2f ar da a dr + f rr dr 2 , (3.2) 

can be decomposed into 5'0(p)-algebraically-irreducible terms in the form 

fab = 0u a u b sx(r) + -9P ab s 2 (r) + a (a M 6) Vi(r) + a ab t(r) , 
p 

far = 9u a s 3 (r) + a a v 2 (r) , (3.3) 



3 It would affect, though, the system in which the wall is dynamical and R varies along the 
worldvolume. 
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where 

9 = V a u\ a a = u b V b u a , cr ab = P a c P b d V (c u d) --P ab (3.4) 

p 

are respectively the expansion, acceleration, and shear of the flow of u a in the metric 
i] ab . Since these terms are algebraically independent, each of the sets of functions, 
Sj (scalar sector), Vj (vector sector), and t (tensor sector) decouple from the others 
in the linearized equations and can be studied separately. 

The Einstein equations R r a = do not involve the / M „ and are independent of 
r. Thus they are 'constraint equations', and can be written in the form 

V a lnr™ +1 = Ua + (n + l)a a . (3.5) 

These equations allow to eliminate the derivatives of r Q (a) in terms of velocity 
gradients. Below we will return to their interpretation in fluid-dynamical terms. 

In order to specify boundary conditions at r = R, we demand that the induced 
metric remains fixed and uncorrected to the order we are working, 

Kb = -UaU b + P ab + 0{d 2 ) . (3.6) 



This requires that 



s 1 (R) = s 2 (R) = Vl (i?) = t(R) = . (3.7) 



In addition, we ask that the stress-energy tensor is in 'Landau frame', defined such 
that the c( 
tions, i.e., 



that the corrections to the leading order value lie entirely along spatial direc- 



u a f^ = 0. (3.8) 

A brief calculation shows that this implies the conditions 

n 4- 1 

vi(i2)=0, s' 2 (i?) = -^s^R) . (3.9) 

The construction of the solution to the Einstein equations for that satisfies 
these boundary conditions and in addition is regular at the horizon, is done in 
appendix |A| using the results of [E]. The explicit results are in eqs. (A. 4). This 



provides the complete metric for the fluctuating black brane, to first-derivative order, 



for any solution of the equations (3.5) 



The solution, however, is written in terms of the velocity field u a and the connec- 
tion V a for the metric r] ab , which is not the physical metric on the wall at r = i?F] 



4 It is neither the metric on the surface at r — > oo, since with our boundary conditions the 
functions Si 2 and t do not vanish there. 
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Nevertheless, we can readily find the relation of the latter to quantities on the wall. 
We do this in appendix [Bj where we find that 

a a ^ a °ab 

U = r^r^r > = r^^r > a ab 



vW yfTWY vW 

a a = a a + -^==P a b d b ^f{R). (3.10) 



The change in the velocity, expansion and shear in (3.10) is simply a local redshift. 
The acceleration is not redshifted, but it is affected by the spatial variation of the 
redshift along the worldvolume. The point is clearer if we introduce the Newtonian 
potential (f), 

f(R) = e 2 *, (3.11) 
which depends on a a through r$. Its spatial gradient is 

V a = P o 6 cW>, (3.12) 

and we see that the modification of the acceleration is due to a 'force' term, 

a a = a a + Va4>. (3.13) 

Now using these relations we write the metric in terms of wall quantities as 

ds2 = (~ UaUb + Pab + J== u a u b Si (r ) + ySS. §p ab S2 ( r ) 

V f{R) v/CR) p 



+ 



^=(a-V<P) {a u b) vi(r) + v / f(R)a ab t(r^jd(T a da b 



VJW) 

+2 (§u a s 3 (r) + (a a - V a 0) v 2 (r)) da a dr 

+ (! + V / 7CR)^ 4 (r)J +r 2 dn (n+1) . (3.14) 

As we will see in the next section, the most consequential effect is the modification 
of the acceleration. 



4 Hydrodynamics and stability 

For later reference we review briefly some generic features of the dynamics of perfect 
fluids. 
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4.1 Perfect fluid dynamics 

The hydrodynamic equations V a T ab = for a generic relativistic perfect fluid with 
stress-energy tensor 

T ab = eu a u b + PP ab (4.1) 

are 

u a u b V b e + P ab V b P +(e + P) {9u a + a a ) = . (4.2) 

The fluid is assumed to satisfy the local thermodynamical relations e + P = Ts and 
de = Tds. Defining also 

v!=§, (4.3) 



we can write the fluid equations (4.2) in a conveniently simple form 



V alias = 9u a - —a a . (4.4) 

Consider now a fluid state initially in static homogeneous equilibrium in its rest 
frame, and introduce a small perturbation, 

s ->• s + 6s e *"*+*-* , u a = (1,0) -)• (l, ( W aJ ' +ikoc ). (4.5) 



Then the solution to the linearized eqs. (4.4) gives fluctuations with dispersion re- 
lation 

u{k) = ^ s k + 0{k 2 ) , (4.6) 

where k — |k|. Hence v s is the velocity of propagation of small density fluctuations, 
i.e., the speed of sound. 

4.2 Black brane hydrodynamics 

In the black brane fluid, the entropy density s is directly related to the horizon 



thickness tq by (2.13). Therefore, density fluctuations in the fluid are variations of 
the horizon radius. 

We can immediately see the hydrodynamic Gregory-Laflamme instability in the 
simplest case in which the cavity wall is removed, R — > oo [U |S]. In this case the 



induced metric is r] ab , the effective fluid velocity is u a and, comparing to (4.4), we 



see that the constraint eqs. (3.5) are the equations of the effective relativistic fluid 



at asymptotic infinity. The effective speed of sound is 

1 



•* = -iTT (47) 
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which is imaginary and therefore fluctuations of Tq grow exponentially in time instead 
of oscillating as sound waves. This is the Gregory-Laflamme instability in the regime 
of long wavelengths and small (imaginary) frequencies. 



Finite cavity: effect of redshift gradients. When we insert the cavity wall 



at finite R, the gradient term in (3.13) modifies the acceleration with which the 
effective fluid responds to a change in r . A local fluctuation 5r > results in a 
smaller 0, which tends to push the effective fluid away from the region of increased 
tq. Conversely, a region of locally smaller gives a gradient term that accelerates 
the fluid towards that region. Therefore, as a consequence of fixing the metric on 
the cavity wall, the creation of inhomogeneities along the worldvolume is hindered. 
The result is to make the fluid more stable. Moreover, the effect is more pronounced 
as the cavity radius R gets closer to the brane, since the redshift becomes stronger]^] 



We can be more quantitative if we use (3.5) to write 

P a b V b \nr = a a (4.8) 

and then 

V »*=-27(ii)-k^V i r„ = --( 7M -l)a a . (4.9) 

Since f(R) < 1, we see that V a is directed opposite to a a and therefore opposes 
the unstable growth of inhomogeneities. If, by decreasing R, the gradient grows to 
a value such that 

V a( f) = -a a , (4.10) 

then in this state the acceleration of the fluid on the wall vanishes, a a = 0: the black 
brane does not react to a density fluctuation, and the instability disappears. This 
happens when 

Km- 1 )^ (411) 

that is, when 



I o\ l/ n 

n + 2 x 



R = R c = r [-^) ■ ( 4 - 12 ) 

If we reduce R below R c , the acceleration a a will be directed against the inhomo- 
geneities, and the black brane will be stable. 

5 Note that the effect is the opposite of what would occur to a material fluid localized on a brane 
at finite R: this would be gravitationally pulled towards larger local mass densities, i.e., larger 
ro. Instead, our effective fluid is not any matter in the spacetime, but rather it is a 'holographic' 
description of the black brane. 
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Effective fluid equations and speed of sound. We can frame this discussion 



in more fluid-dynamical terms. From (4.9), the relation between the accelerations 



(3.13) in the black brane fluid is 



a a = a a V \, (4.13) 



where 



v 2 



n ( 1 
n + 1 V~2\J(R) 

1 1 - (Rc/R) n 
n + 1 f(R) 



(4.14) 



Now the constraint equations (3.5) written in terms of the effective fluid velocity 
on the wall become 



1 



W a \nr^ L =6u a --a a . (4.15) 

vi 



Comparing to the general form of the perfect fluid equations (4.4), we see that (4.15) 



are the equations V a T = for the stress-energy tensor (2.9) and v s = {dP/de)R 



is indeed the speed of sound. The quasilocal stress-energy tensor is known to be 
conserved on general grounds |17] . so the result is not surprising. What we have 
done here is to see explicitly how these conservation equations emerge, on walls at 
finite R, from the Einstein constraint equations. 



The explicit form of eqs. (4.15) is in any case very illustrative. They show clearly 
that, to the order we work, all the flow with R of the black brane dynamics is due 
to the modified acceleration term. Since the entropy density is independent of R, 
the change of the effective fluid with R can be fully accounted for by the change of 
v 2 



s ' 



The sound velocity v s is imaginary for large R, but it vanishes when the cavity 



reaches the critical radius R c in (4.12), and then becomes real for cavity radii R G 
(tq, R c ). This change in stability works in the direction expected from our argument 
above. 

Intuitively, the speed of sound is a measure of the rigidity of the system to 
worldvolume fluctuations. For very large cavities, the geometry is excessively soft, 
indeed 'anti-rigid', to the point of being unstable to deformations. The cavity wall, 
by fixing the geometry at a finite distance from the black brane, increases its stiffness 
and can even render it stable when the cavity is small enough. In fact, v 2 grows 
without bound as the wall approaches the horizon, thus making the effective fluid 
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incompressible in that limit J^] 

The growing stiffness of the system caused by worldvolume gradients of the 



redshift is also apparent in the expressions for the extrinsic curvature, (2.8), and 



the pressure, (2.12). Thus, the effective hydrodynamic theory explains in a simple 



manner why and how the black brane turns from unstably soft to stably stiff. 



Unstable perturbation. It is now easy to give the complete form of the unstable 
black brane solution in the cavity. We illustrate it, for simplicity, in the case of a 
black string with worldsheet coordinates a a = (t,x). Take a velocity profile of the 
form 

-1, u x = exp v l ktj cos(kx) 5u , (4.16) 



and work to linear order in the small amplitude Su. The metric is given by (3.14), 



with the functions Sj(r), v.;(r) as in eqs. (A.4), and with 



(4.17) 



9 = —k exp C\/—v% kt \ sin(kx) 5u , 
ix = V~ v s kexp (^\/—Vg kt^j cos(kx) 5u . 



and 



f(R) 



l-(R c /R)n X - 



(4.18) 
(4.19) 

(4.20) 



4.3 Correlated dynamical and thermodynamic stability 

We can easily see that the change in dynamical stability at R = R c corresponds 
precisely to the change in the local thermodynamic stability of the black brane, i.e., 
in the thermodynamic stability of the black hole that one obtains at any given point 
on the worldvolume. Since (2.15) and (2.16) imply dP = sdl~ we have 




s 



(4.21) 



where Cy is the specific heat at fixed volume. Since the sign of Cy determines the 
local thermodynamic stability, the connection between the latter and the dynamical 



6 The fact that v s — > oo as R — > tq does not necessarily entail any violation of causality. 
Hydrodynamic fluctuations are low-frequency modes, and causality is controlled by modes in the 
high-frequency end of the spectrum, see e.g., the discussion in [15] . 
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stability of the brane, in the hydrodynamic regime, is obvious. This is nothing but 
the fact that stability of hydrodynamic modes associated to conserved quantities is 
governed by the local thermodynamic properties of the fluid. 



In the calculation of the speed of sound and in eq. (4.21) we only need the static 



brane solution. What our study of the fluctuating brane shows is that there is indeed 
an explicit solution for a black brane in a fixed cavity which is regular on the horizon 
and which is dynamically stable or unstable in accord with its thermodynamical 
stability. 

5 Viscous hydrodynamics 

Having the fluctuating black brane geometry to first order in velocity gradients, we 
extract its quasilocal stress-energy tensor at finite R including dissipative terms. 



5.1 Bulk and shear viscosities do not run 

The general form of the stress-energy tensor on the cavity wall at finite R, in the 



spacetime given by (3.1), (3.2), (3.3), and with boundary conditions (3.7) and (3.9), 
is 



T, 



ab 



[n 



8n~G L 
-C9P ab -2fja ab + 0(d 2 



l)R n ^J(R)u a u b + d R (R n+1 ^/J{Kj) P t 



ab 



with bulk and shear viscosities 



fin 



+ 1 



8ttG 



R 



n+l 



<(R) 



and 



O tfn+1 



(5.1) 



(5.2) 



(5.3) 



8ttG 4 

Substituting the explicit values for the solution that is regular on the horizon, we 
get 

e /1 1 \ 

(5.4) 



c = ^(i + ' 



and 



2tt \p n+l 



An 



(5.5) 



with s the entropy density (2.13) 
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While the result for the shear viscosity is not surprising, the fact that the bulk 
viscosity remains the same at all R is probably less obviously expected. In particular, 
observe that it is of the form 

< = £ G - • (6 - 6) 

Thus it depends on the asymptotic value of the speed of sound, instead of its value 
v 2 at the cavity wall, which one might naively have guessed. Had it been the latter 
case, ( would have run with R. Since s is independent of R, we can equivalently say 
that neither fj/s nor £/s run with R. 

This absence of running of £ is most probably related to the one in [16] . where it 
was found that for AdS black branes in a finite cavity the bulk viscosity remains zero 
at all R, despite the fact that the wall breaks conformal invariance. There exists an 
explicit mapping between AdS gravity and the sector of vacuum gravity involved in 
our system [19] which is independent of the cavity wall. Conceivably, it relates our 
result to that of [TB] and possibly makes clearer why ( depends on v 2 instead of v 2 s . 

For AdS black branes, ref. [TB] found an intriguing relation, T a a = —de/dlnR, 
for the running of the energy density with R. For our neutral black branes, the 
same equation formally applies if we set n = —p — 1, which is not any physical 
black brane, and in fact corresponds to setting D = 2. Since, again formally, when 
n = — p — 1 one gets £ = 0, this running is valid including first-derivative corrections. 
The reason why this result holds in this context is possibly related to properties of 
analytic continuation in n, but at present its ultimate meaning is unclear to us. 



5.2 Spectrum of unstable modes 



With £ and fj we can compute the corrections to the fluid equations due to the 
viscous damping of density fluctuations. This gives us a better approximation for 
the spectrum of unstable modes at finite R. 

Solving the fluid equations to quadratic order in momenta k, the unstable modes 
of the black brane in a cavity with R > R c have imaginary frequency 



fi(Jfe) = ^/^k- 



2Ts 



1 277 



c 



k 2 + 0{k 3 



(R c /R) r 



k 2 r - 



n + 2 



VJiR) + o{k 2 



(5.7) 



v^TTV f{R) "n{n + l) 

where R c is given in (4.12). We illustrate this result in fig. [I] As R approaches R c , 
the instability gets weaker, having both a smaller rate of growth Q and a shorter 
range of unstable wavenumbers k. 
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kr 



Figure 1: Spectrum of GL unstable modes for a black brane in a cavity of radius R, 



according to eq. (5.7). The curves correspond to n = 1, for which the critical radius is 
R c = 1.5 r . 



While there is no previous calculation of black brane instabilities in a finite 
cavity that we can match these curves to, we can compare against the computation 
in ref. [13] of the wavenumber kci of the zero-mode, for which Q(kcL) = 0. Our 



analytic expression (5.7), truncated to quadratic order, gives 



1 ny/n+~T y/1 - (R c /R) n 
n + 2 



(5.* 



r n + 2 f(R) 

We display this result in fig. [2j where we compare it with the corresponding one 
of Fig. 3 in [13]. The qualitative agreement between the two graphs is apparent, 
but one can easily discern quantitative discrepancies. These are expected, since 
( 5.8[ ) has been obtained under the hydrodynamic assumption of small wavenumbers 
k/T <C 1, which is not satisfied in general. As in [8], we may expect the agreement 
to improve for larger n. 



5.3 Critical behavior 



Eq. (4.15) implies that when v s = 0, the acceleration of the fluid under a density 



perturbation vanishes, i.e., the fluid does not respond to variations of tq. Using 



(4.21) we see that at that point the effective specific heat Cy becomes infinite. 



This result makes manifest that the divergence of the specific heat is linked to the 



16 



<k GL r ) 2 




12- 



10- 



Pb o 




Figure 2: Evolution of the GL zero mode koL as a function of the cavity radius R. Left: 



obtained from the approximate analytical expression (5.8). Right: numerical results from 
[15] . Following [15] . in the vertical axis we display the 'tachyon mass squared', —{kcLfo) 2 , 
and in the horizontal axis pf, = (R — ro)/{R c — vq), so that the critical radius is always at 
Pb = 1. The curves are for n = 2, . . . , 7, and lie closer to the horizontal axis the smaller n 
is. The finite, non-zero slope of the curves at pb = 1 indicates mean-field critical behavior 
(Pi- 



ghost instability — the hydrodynamic mode is at the threshold of becoming ghost- 
like. The connection between the divergence of Cy and the tachyonic instability 
is, instead, only indirect: as discussed in the introduction, a hydrodynamic ghost 
instability of a black hole system is typically accompanied by a tachyonic instability. 
The hydrodynamic statement that dP/di — > means that the fluid does not react 
with any pressure gradient to local variations of the energy density. Thermodynam- 
ically, the fact that Cy — > oo means that under these density variations the system 
does not create any temperature gradients that would restore it back to thermal 
equilibrium. Both effects, hydrodynamic and thermodynamic, are of course related 
via dP = sdi~. At the point where Cy — > oo, thermal fluctuations in the fluid have 
infinite wavelength. Then, this is a thermodynamic critical point. 
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In the critical state, the tachyonic zero mode becomes massless, i.e., has infinite 
wavelength, and therefore it must be included alongside with the hydrodynamic 
modes in the effective low-frequency, long wavelength theory. Note, however, that 
this tachyon is never a proper hydrodynamic mode: the fluid equations never admit 
(to any arbitrary derivative order) a non-trivial static solution. The inclusion of the 
massless tachyon in the effective theory at the critical point cannot be done working 
solely within hydrodynamics. 



Thus, the result (5.8) from a truncated hydrodynamical calculation need not be 



accurate near the critical point. Note, however, that it predicts a critical behavior 

k GL ~{R-R c ) l '\ (5.9) 

which has a mean-field theory exponent. This appears to be correct: the numerical 
curves of [13] for — fc^ L , reproduced in the right-side graph in fig. [2l cut the horizontal 



axis with non-zero, finite slope consistently with the critical behavior (5.9). Perhaps 
we should not be too surprised: a mean-field theory description of the critical state 
might be natural in a classical field theory like general relativity and in a state with 
smooth horizons. 

First order transition: inaccessibility of the critical point. The locally- 
unstable phases of the system of a black brane in a cavity are extremely fine-tuned 
states, which we have considered not so much for their possible relevance to ac- 
tual physical phenomena but as being illustrative of the dynamics of black branes. 
But, actually, even the critical point and the associated second-order phase tran- 
sition cannot be expected to be reached through any physical process, as they are 
dominated by a stronger first-order transition. 

We have presented the problem starting from the black brane inside a very large 
cavity and then reduced the size of the cavity until the black brane becomes locally 
stable. However, from the perspective of a physical process it makes more sense 
to start from the system in the stable regime of R <C R c and then follow it as we 
increase R. We are imagining that we keep the horizon size r fixed, but one could 
easily consider other processes, e.g., fix R and change the temperature of the box. 

When R reaches the value 

/ n + 2 \ 1/n 

fll =U^T)) W ' (510) 

the pressure P of the black brane becomes equal to that of Minkowski space in 



the same cavity, (2.19), and for R > R\ the latter has larger pressure. Thus, 



18 



when the brane is 'supercooled' at R± < R < R c , even if it is locally stable it 
will undergo a first-order phase transition and spontaneously nucleate bubbles of 
the 'true vacuum', i.e., hot flat space^J This will form holes in the worldvolume of 
the brane that then begin to expand. It may be interesting to study further this 
phenomenon, which shares features with the one studied in [20]. At any rate, it 
will prevent the observation of the critical state at R = R c , and indeed of all the 
unstable black brane phases. 
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A Solution with finite-cavity boundary conditions 

Ref. [S] solved the Einstein equations for the perturbations to first order in deriva- 
tives and obtained the general solution that is regular on the future event horizon. 
The solution, expressed in Eddington-Finkelstein coordinates, contains a number of 
integration constants that are to be determined by conditions at the spatial bound- 
ary. For ease of comparison, we use the same notation for these constants (c vr , Ca, 
etc.) as in [S|. All the calculations at finite R, including the stress-energy tensor, 
could be carried out in Eddington-Finkelstein coordinates [6], but here we choose 
to work in Schwarzschild coordinates to maintain continuity with [5], [8] . In these 
coordinates the solution is 

s i( r ) = o/ W + \ \f( r ) ) s 2( r ) + °w - 2c w - f{r)c vv , (A. la) 



2(n + l) V n + 2 

s 2 (r) = ^ln /(r)+c«, (A.lb) 
n 



7 For the black hole inside a spherical cavity (p — 0) this is the analogue of the Hawking-Page 
transition. 
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2/(r) 



" /(r) ) s 2 (r) + (/(r) 



n + 1 

-2c vr 



n + 1 



( 




n 




V r o 





f(r) 



s 4 (rj 

vi(r) 
v 2 (r) 



, ^ r (1) 
r* — r + vi(r) 

S2(r) + Cjj , 



2(n+l) 



where 



dr 



(A.lc) 
(A. Id) 

(A.le) 
(A.lf) 
(A.lg) 

(A.2) 



The function f r j(r) is a gauge-dependent function only constrained to be finite on 
the horizon and thus could be set to zero. 



The boundary conditions (3.7), (3.9) are satisfied by choosing 

2r 



(■-il. 



n 



\nf(R) 



with which 
Si(r 
s 2 (r 

s 4 (r 

vi(r 
v 2 (r 

t(r 



n + 2 
2(n + 1) 
2r f(r) 



fir) 



n 



n + 2 



s 2 (r 



n f(R) 



n 



2/(r] 



2r n 



, -/(r) U(r) + — ^-(/(r)-l) n- + 1 
n + 1 / n+1 V^o 



4r - ns 2 (r] 
2(n + 1) 



0, 

— r 
s 2 fr) . 



+ /«(r) 



(A.3) 

(A.4a) 
(A.4b) 
(A.4c) 

(A.4d) 

(A.4e) 
(A.4f) 

(A.4g) 



This is the unique solution, up to the gauge choice of f r j(r), that satisfies the 
regularity condition at the horizon and the boundary conditions on the wall at fixed 
R. 
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B Connection on the wall 

We present here the relation between the connections V a and V a compatible with, 
resp., the metrics r\ ab and h a b, and we use this to relate the gradients of the respective 
velocity vectors, u a and u a . We decompose these gradients, as usual, into traceless 
symmetric shear a ab , expansion 8, acceleration a a , and antisymmetric vorticity u ab , 
so that 

V a M 6 = (Tab + -OPab ~ U a a b + U ab , (B.l) 

V 

and similarly for hatted quantities. 

We follow the same steps as in [16]. If the difference between the metrics is 



lab = Kb ~ Vab (B.2) 

then the difference between the connections, T ab , such that 

V a V b = V a V h - f c ab V c (B.3) 

is given by 



Kb = ^ {^albd + Vblad ~ Vdlab) ■ (B.4) 



In our case 



lab = (1 - /) u a u b . (B.5) 

Here we always take / evaluated on the wall, i.e., f = f{R), which depends on a 
through r (cr). We find 



-, c u d d d f 1 d 1 

lb = -^pu c u a u b - jU c u (a P b) d d d f + -1 



+ - 1) u c L ab + UPabj + 2(1 - f )u (a uj b) c - (1 - f)a c u a u b . (B.6) 



Using this and (2.6) we find that 



V a w& = (°ab + + \f]u ab - u a (^a b + ^P b dJ ) , (B.7) 



from where we immediately deduce eqs. (3.10) and 



&ab = Vf^ab- (B.8) 

Note that we have not made use anywhere of the fluid equations of motion for 
eliminating the derivatives of / in favor of derivatives of the velocity. This is done 



in sec. 14.21 
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An alternative but equivalent way of obtaining the same results is the following. 
We may regard the two metrics as related, to zeroth derivative order, by the change 
u a a a = a/ f(R) u a a a . Above we have set u a = y f(R) u a and left the coordinates 
unchanged. But we could just as well leave u a = u a , perform a coordinate rescaling 
(of time), and include the derivatives of f{R) that result from this coordinate change 
into the correction terms f^ u . 
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